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Summary 
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1.  Potential  flow  theory  predicts  that  a body  accelerated  or 
decelerated  in  motion  in  a straight  line  in  a fluid  will  experience 
an  inertia  force  greater  than  thet  due  to  the  mass  of  the  body  - a 
so  called -%dded  mass y effect.  Submarine  trials  confirm  the 
phenomenon  but^indicate  tfyat  for  motion  in  the  fore  and  aft 
direction  the  Padded  massJ  is  much  greater  than  predicted.  This 
could  be  due  to:-s 


J 


Cx-’tyCt)  Transient  loss  of  thrust  from  the  propeller^ 

- — Changes  in  skin  friction  due  to  the  acceleration. 

J 

^>(i)^£hanges  in  form  drag  due  to  the  acceleration.  _ 

2.  lowing  tests. have  been  conducted  on  a body  without  propeller. 
These  eliminate  s(‘  as  tbe  primary  cause.  - 

^ ^ & 

3.  A theoretical  and  experimental  study  of  has  been  made  for 

a flat  plate.  Theory  predicts  an  increase  in  skin  friction  for 
both  laminar  and  turbulent  boundary  layers,  but  experiment  on  the 
‘turbulent  layer  only  partially  supports  the  theory.  Although  further 
analysis  of  the  data  is  desirable,  the  results  suggest  that 
acceleration  may  restore  a measure  of  laminar  flow  in  the  boundary 
layer.  Further  work  should  concentrate  on  this  aspect  and  on  the 
effect  of  acceleration  on  form  drag. 


Approved  for  issue 
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THE  DRAG  OF  MY  ACCELERATING  SUBMARINE 
PART  I SKIN  FRICTION 


By  T B B.oth 
1.  INTRODUCTION 

An  accurate  knowledge  of  the  acceleration  and  deceleration  of  a 
submarine  has  become  an  important  aspect  of  the  calculation  of 
safety  trajectories  - a subject  thrown  into  prominence  by  the 
advent  of  the  high  speed  submarine.  Without  this  knowledge  it  is 
not  possible  to  calculate  accurately  the  depth  from  which  a 
submarine  can  reach  the  surface  in  the  event  of  flooding,  neither  is 
it  possible  to  calculate  accurately  the  depth  excursion  following  a 
hydroplane  jam. 

It  is  well  known  that  a body  accelerating  in  a fluid  appears  to 
experience  an  inertia  force  greater  than  (mass  of  the  body)  x 
acceleration.  This  apparent  inertia  force  can  be  expressed  in  the 
form  ('virtual  mass')  x acceleration,  the  difference  between  the 
acual  mass  and  tho  virtual  mass  being  commonly  referred  to  as 
'added  mass'.  The  terms  'virtual  mass'  and  'added  mass'  are 
unfortunate,  being  historically  associated  with  the  concept  of 
'entrained  water'  ie  water  dragged  along  with  the  submarine  (in 
addition  to  free  flooding  water  trapped  within  the  form  which  is 
legitimately  part  of  the  dynamical  mass  of  the  submarine).  The 
increase  in  inertia  is  in  fact  due  to  the  need  to  increase  the 
kinetic  energy  of  the  fluid  whenever  the  velocity  of  the  body  is 
increased.  Lamb  (Reference  l)  obtained  theoretical  solutions,  based 
on  the  work  of  Green  (Reference  2),  for  the  kinetic  energy  stored  in 
the  ir.viscid  potential  flow  round  an  ellipsoid.  He  showed  that  the 
excess  of  the  inertia  force  over  (mass  of  the  body)  x acceleration, 
could,  for  acceleration  along  a principal  axis,  be  expressed  in  the 
form 


k x (mass  cf  fluid  displaced)  x acceleration 

where  k is  a constant  whose  value  depends  on  the  axis  concerned, 
and  which  he  refers  t.r>  as  an  'inertia  coefficient'.  He  obtained 
similar  coefficients  for  angular  accelerations.  The  'added  mass' 
and  the  'inertia  coefficient'  are  related  by: 

added  mass  = k x (mass  of  fluid  displaced) 

Note  that  the  inertia  coefficient  factors  the  mass  of  fluid  displaced 
ie  the  form  displacement,  which  is  not  necessarily  the  same  as  the 
mass  of  the  body.  Reference  3 gives  the  values  of  the  inertia 
coefficients,  calculated  according  to  Lamb's  theory,  for  a wide  range 
of  proportions  of  ellipsoids. 

The  limitations  of  potential  flow  theory  are  well  known,  and  until 
the  invention  of  the  Planar  Motion  Mechanism  (References  ^ and  5) 
there  was  little  experimental  data  concerning  the  inertia  coefficients 
of  a body  in  water.  We  now  know  that  Lamb's  predictions  of  the 
inertia  coefficients  for  vertical  and  lateral  motions  are  of  the 
correct  order  when  applied  to  modern  submarine  shapes,  and  may  be  as  high 
as  0.9.  There  is  no  doubt  that  for  these  motions,  inviscid  potential 
flow  accounts  for  nearly  all  of  the  increased  inertia. 


1 


When  applied  to  the  forward  motion  of  a submarine,  Lamb's  theory 
predicts  very  small  values  of  k,  typically  less  than  0.05  and  this 
as  we  shall  see  is  far  too  low.  It  is  this  discrepancy  which  is 
the  subject  of  this  note. 


The  first  computations  of  submarine  acceleration  and  deceleration  at 
AFW  were  carried  out  some  ten  years  ago,  and,  in  common  with 
practice  elsewhere,  the  inertia  coefficient  for  fore  and  aft  motion 
was  based  on  that  predicted  by  Lamb,  replacing  the  submarine  by  an 
equivalent  spheroid,  while  propeller  thrust  was  based  on  quasi-steady 
state  values.  (There  is,  in  fact,  sor.e  justification  for  this 
practice  in  Reference  6,  in  which  the  inertia  coefficient  was  obtained 
experimentally  from  a drop  test  in  water  of  a body  of  fineness 
ratio  3.5.  The  v&lue  of  k was  deduced  to  be  O.OL,  which  is  less  than 
the  theoretical  value  based  on  Lamb.) 

Some  full  scale  data  became  available  in  late  1967  in  the  form  of 
acceleration  and  deceleration  runs  of  HMS  VALIANT  for  which  the 
velocity  and  rpm  were  recorded  against  time.  Although  not  entirely 
satisfactory,  and  having  rpm  time  histories  wh^ch  differed  from  those 
assumed  in  the  predictions,  it  was  evident  that  the  submarine  neither 
accelerated  nor  decelerated  as  quickly  as  predicted. 


Opinion  at  that  time  attributed  the  discrepancy  to  the  transient 
conditions  in  which  the  propeller  was  working  compared  with  the  quasi- 
steady condition  assumed  in  the  prediction.  Nevertheless  the  runs 
were  analysed  assuming  that  quasi-steady  conditions  did  apply.  Thrust 
and  resistance  were  estimated  using  the  steady  state  propeller 
characteristic  and  ITTC  skin  friction  values,  together  with  values  of 
'Taylor  Wake',  'Augment'  and  speed  rpm  ratios  obtained  from  speed 
trials,  and  from  them  the  'virtual  mass'  was  deduced.  This  work  is 
reported  in  Reference  7 from  which  Figure  1 is  taken. 


VIRTUAL  MASS 

FORM  DISPLACEMENT 

1-5*. 


♦ 


• accelg 

♦ decelg 


- ♦ 


♦ 


JL. 

10 


J ■ — . > 


speed  kts 


-I 

20 


Figure  1 . Trials  Result 


There  is  some  scatter  on  the  results,  but,  assuming  the  submarine  to 
have  been  approximately  neutrally  buoyant,  an  inertia  coefficient  of 
between  0.15  and  O.h  is  indicated.  This  is  too  high  to  be  accounted 
for  by  the  error  incurred  by  treating  the  submarine  as  an  equivalent 
ellipsoid. 

The  high  value  of  inertia  coefficient  could  be  due  any  of  the 
following: 

1.  Transient  conditions  on  the  propeller  and  associated 
effects  on  the  hull  causing  a loss  of  thrust. 

2.  Effects  of  acceleration  and  deceleration  on  the  skin 
friction  of  the  hull. 

3.  Effects  of  acceleration  and  deceleration  on  the  form 
drag  of  the  hull. 

1 . 1 . Propeller  Transients 

To  test  the  first  of  these  possibilities  a submarine  model  without 
propeller  was  towed  in  No  2 Ship  Tank.  The  resistance,  X,  was  measured 
during  the  accelerated  and  decelerated  parts  of  the  tow,  as  well 
as  at  steady  speed.  Figure  2.  Subtracting  the  steady  state  drag, 

Dy,  corresponding  to  the  speed  at  any  instant, from  the  total 

resistance  gives  the  component  due  to  acceleration.  Dividing  by 
the  acceleration  at  that  instant  gives  the  'virtual  mass'  of  the 
model.  Figure  2 shows  this  to  be  10  to  30  per  cent  of  the  form 
displacement.  This  result,  typical  of  that  obtained  in  several 
similar  experiments  on  different  models,  indicates  that  the 
propeller  is  not  the  primary  cause  of  the  phenomenon. 

More  recently  the  transient  forces  on  the  propeller  have  been 
measured  (Reference  8)  and  have  been  found  to  be  insignificant.* 


* This  conclusion  is  confirmed  by  early  experiments  on  airship  models, 
References  9 and  10.  The  experimenters  did  not  have  towing  tank 
facilities  and  their  results  are  unsatisfactory.  Nevertheless  they 
confirm  the  increase  in  drag.  The  reason  for  the  increase  appears  to 
have  attracted  little  attention.  Presumably  interest  in  airships  was 
waning  in  favour  of  fixed  wing  aircraft. 
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1.2. 


Form  Drag  and  fkin  Friction 


The  effects  of  acceleration  on  both  form  drag  and  kin  friction  merit 
: avert i nation.  The  latter  has  teen  studied  first,  partly  because  there 
are  rood  rea..ons  to  expect  that  acceleration  will  increase  the  skin 
friction  and  partly  U cause  it  is  more  readily  investigated.  The 
results  are  reported  here. 

The  effects  of  acceleration  on  form  drag  are  less  clear,  and  are  the 
subject  of  future  work,  which  will  in  due  course  form  a Part  II  to 
follow  this  report. 

2.  SKI U FRICTION  oP  A!.’  ACCELERATING  PLATE 

The  derivation  of  equations  of  steady  state  boundary  layer  flow  used 
in  this  section  is  given  in  Reference  11  or  more  conveniently  in 
Reference  12. 

There  is  no  doubt  that  the  skin  friction  of  a flat  plate  is  increased 
above  the  steady  state  value  by  acceleration.  A simplified  physical 
explanation  is  as  follows. 

Suppose  a flat  plate  of  infinite  span  to  be  in  stead'  motion  in  its 
own  plane  and  at  right  angles  to  its  leading  edge  so  that  the  velocity 
profile  in  the  boundary  layer  at  some  point  P is  as  shown  by  curve  A 
on  Figure  3* 


DIR  ECTION 

OF  MOTION  OF  PLATE 


Figure  3.  Velocity  Profile  in  the  Boundary  Layer 


Now  suppose  the  plate  to  be  accelerated  sharply  in  the  same  direction 
An  instant  lat'-r  the  velocity  Las  increase  i ani  the  velocity  profile 
which  would  develop  at  F if  this  increased  velocity  were  held  steady 
is  shown  by  curve  E.  but  tnat  part  of  the  boundary  layer  away  from 
the  surface  of  the  plate  does,  not  have  time  to  develop,  so  the 
velocity  tnere  remains  close  to  curve  A.  At  the  surface  of  the  plate 
on  the  other  hand,  the  velocity  in  the  fluid  must  increase,  on  the 
assumption  that  the  hypothesis  of  'no  slip'  remains  true.  The  actual 
velocity  profile  is  therefore  of  the  form  of  curve  C.  The  velocity 
Gradient  at  the  surface  is  creator  than  in  steady  flow  and  hence  the 
skin  friction  is  creator. 

’ . 1 . Hot  ut i on 


y 


u 


. L_u 

0 77  1 H / 77  YJ'nTTPTTT 


X 


Oxy 


1 b t 
P M V 


u u 


u v 


6 


Figure  b 

Axes  fixed  in  the  plate. 

O origin  at  the  leadinr  «*dce. 
x co-ordinate  along  the  plate, 
y co-ordinate  normal  to  the  plate. 

Lengthj  breadth;  thickness  of  the  plate. 

Density , viscosity * kinoTn?tic  -r  of'  th^ 

fluid . 


Velocity;  acceleration  of  the  undisturbed  stream 
relative  to  the  plate. 

Velocity  components  in  the  boundary  layer  at 
(x  , y)  relative  to  the  plate. 


Boundary  layer  thickness. 


Displacement  thickness 
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L’rug ; Jrac  a’  *onstant  sp<cd  U;  drag  due  to 
acceleration  = D - 1^. 

Local  skin  friction  per  unit  area;  etc. 

Mkin  friction  of  one  side  of  a flat  plate;  etc. 

Mean  skin  friction  coefficient  of  a flat  plate 

(one  side)  = F,,/jpU2  lb. 
u 


3.  LAMINAR  BOUNDARY  LAYER 


3.1.  Equations  of  Motion 


The  equations  of  motion  are  conveniently  established  by  considering 
the  plate  to  be  stationary  in  an  accelerating  stream.  In  order  to 
. accelerate  the  stream  it  is  necessary  to  introduce .a  pressure 
gradient  in  the  direction  of  flow  of  magnitude  - pU  where  U is  the 
acceleration  of  the  undisturbed  stream.  Hence  for  two  dimensional 


(D 


(2) 

The  equations  of  continuity  and  laminar  skin  grief-inn  remain 
unchanged 


motion  the  usual  equation  of  motion 


3u  _3u  _3u  1 j,  32u 

3t  U 3x  V 3y  p 3x  V 2 


becomes 


3u  3u  , 3u  ' . 32u 

— + u — + v—  =U  + v 

3t  3x  3y  2 


3u 

3x 


+ 


(•;) 


('«) 


3.2.  Exact  Solution 


The  equations  have  been  solved  exactly  for  a plate  of  length  It.  6 m 
accelerated  from  rest  in  water  at  0.6  m/s2,  by  the  computation 
described  in  Appendix  I.  After  the  plate  has  travelled  a distance 
of  three  lengths  the  speed  is  held  constant  ( U . 1 m/s).  The  growth 
of  skin  friction  at  three  points,  distant  1/5,  21/5  and  1*1/5  from 
the  leading  edge  is  shown  in  Figure  5.  The  points  at  1/5  and  1*1/5 
were  selected  to  be  free  of  any  error  which  might  arise  nearer  the 
leading  and  trailing  edges  due  to  the  computing  singularities  which 
occur  there. 
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iffure  5.  Skin  Friction  of  the  Laminar  Boundary 


3.3.  Quasi  Steady  Solution 


The  simplest  possible  solution  is  to  assume  that  the  skin  friction  at 
any  instant  is  that  corresponding  to  fully  developed  steady  flow  at 
the  instantaneous  speed.  This  anounts  to  ignoring  the  effects  of 
acceleration  and  cannot  be  expected  to  give  the  correct  value,  but  it 
is  of  interest  since  comparison  with  the  exact  solution  isolates  the 
effects  of  acceleration. 

The  skin  friction  is  given  by  the  well  known  solution  for  steady 
laminar  boundary  layer  flow  (References  1?  and  11). 

ty  = 0.3'32eU2(v/Ux)*  (5) 


The  value  at  the  point  21/5  is  given  in  Figure  5 and  shows  that,  as 
expected,  the  quasi  steady  value  is  an  underestimate.  The  effect  of 
acceleration  dominates  in  the  early  stages,  reducing  as  the 
acceleration  continues,  the  effect  being  greatest  on  the  rear  parts 
of  the  plate.  It  is  clear  that  the  increase  of  skin  friction  due  to 
acceleration  is  not  constant  and  that  the  concept  of  a constant 
'added  mass'  does  not  apply  to  viscous  effects. 

3.^ • Blasius  Solution 

* 

An  approximate  solution  for  constant  acceleration  from  rest  was  given 
by  Blasius  based  on  the  omission  <^f  the  convection  term  u from  the 

dX 

equation  of  motion  (1). 

+ (2n2  + l)erf  n)  (6) 


da  (T)  (8) 


the  distance  from  the  leading  edge 
skin  friction  given  by 

2U(py t/v)i  (9) 

\ - /y  = 0 

The  Blasius  solution  agrees  very  well  with  the  exact  solution  in  the 
early  stages  of  the  motion  (Figure  5)-  As  the  acceleration  continue;, 
agreement  remains  good  at  the  point  x = 41/5,  but  underestimates  skin 
friction  at  tne  forward  points,  being  about  !)0  per  cent  low  at  the 
point  x = 1/5  at  tne  end  of  the  accel  rated  motion.  Clearly  the 
convection  terms  omitted  in  the  Blasius  approximation  can  have  an 
important  effect  on  the  forward  part  of  the  plate. 

A physical  explanation  of  the  discrepancy  is  as  follows.  The 
Blasius  solution  allows  for  the  development  of  the  boundary  layer 
by  diffusion  only,  leading  to  uniform  conditions  over  the  plate. 

But  the  boundary  layer  also  develops  by  convection  from  the  leading 


u = U {-  2n2  + 2ti  2-  n 


1 ? 

— n / 

n e 


where  n = (y/2)/(vt)2  erfn  = 


The  solution  is  independent  of  x, 
of  the  plate  and  leads  to  uniform 


( 3u 

= 


= 2pU//7TVt  = 


edge  where  the  boundary  layer  is  thin  and  the  skin  friction  high. 

The  kinetic  energy  required  for  convective  development  is  provided 
by  the  work  done  by  increased  skin  friction.  However  in  the  initial 
stages  of  the  accelerated  motion  the  effects  convected  from  the 
leading  edge  do  not  reach  points  downstream  until  those  points  have 
travelled  approximately  to  the  initial  position  of  the  leading  edge. 
Only  up  to  this  time  can  the  Blasius  solution  be  expected  to  apply. 
The  times  to  travel  l/5>  21/5  and  Ul/5  are  marked  on  Figure  5 and 
support  this  explanation. 

3.5*  Momentum  Theory 


An  alternative  solution  may  be  obtained  from  the  momentum  equation 
x/pU2  = ~ + (1/U)  (00  + A 1 ) + (1/U2)  ~ (U6 1 ) (10) 


where  ^ = 


(1  - u/U)dy 


(11) 


0=f  (1  - u/U)  (u/U)dy 

‘ ) o 


(12) 


3U 

and  in  which  the  terms  in  — vanish  in  this  application. 


In  steady  state  flow  the  following  relations  are  obtained 
<$,  = 1.72  x(R xf  2 


(13) 


0.332  R 


where  R = Ux/v 
x 


(Hi) 


and  we  may  assume  that  these  relations  still  hold  when 
accelerated.  On  this  assumption,  for  the  accelerating 

x/pU2  = 0.332  R “ 5 + (1/U2)  ~r  (1.72  x (v/Ux)sU) 


= 0.332  Rx“  5 {1  + 2.59  (Ux/U2 ) } 


the  plate  is 

pi 


(15) 


Note  that  it  has  not  been  necessary  to  assume  that  the  acceleration 
is  constant. 


The  first  term  is  that  due  to  quasi  steady  flow  while  the  second  is 
due  purely  to  the  acceleration.  The  contribution  of  the  latter  to 
t is 

1 1 , 

t.  = 0.332  R ? (2.59  (Ux/U2 ) }pU2  = O.BGpx  R " "U  (1(>) 

U x x 

Comparison  with  the  exact  solution  in  Figure  6 
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shows  remarkably  good  agreement  in  the  later  stages  of  the  motion. 

It  appears  that  agreement  becomes  good  at  a particular  point  along 
the  plate  when  that  point  passes  the  initial  position  of  the  leading 
edge.  Before  this  occurs  the  momentum  theory  overestimates  the  skin 
friction. 


This  is  consistent  with  the  discussion  of  the  Blasius  solution. 

The  values  of  6 1 and  t used  in  conjunction  with  the  momentum 

equation  are  essentially  based  on  a fully  convected  flow,  leading  us 
to  expect  the  agreement  which  is  in  fact  obtained. 

4.  TURBULENT  BOUNDARY  LAYER 

Despite  the  considerable  volume  of  literature  on  turbulent  flow  the 
author  is  not  aware  of  any  simple  'law  of  turbulence'  which  can  be 
used  to  calculate  the  effect  of  acceleration  on  skin  friction. 
However  the  results  of  applying  the  momentum  equation  to  the  laminar 
boundary  layer  suggest  that  a similar  theoretical  approach  might 
be  made  to  the  turbulent  boundary  layer.  As  a test  of  the  results 
it  is  possible,  although  not  easy,  to  measure  the  skin  friction 
experimentally  on  a plate  towed  at  accelerating  speed  in  the  towing 
tank. 


4.1.  Momentum  Theory 

Equation  (10)  is  valid  for  both  laminar  and  turbulent  boundary  layers 


T/pu2  = If + (1/u)  If  (2Q  + V + (1/u2)  ft  (UV  (18) 


in  which  as  before 


au 

3x 


= 0. 


For  Reynolds  Numbers  of  the  order  of  106  the  'seventh  power  law' 


u/U  = (y/6) 1/7 

(19) 

is  known  to  give  a good  approximation  to  the 
boundary  layer  leading  to  expressions  for  6^ 

_ 1/c 

6 = 0.0463  xRx  /5 

velocity  in  the  turbulent 
and  0 (Reference  12) 

(20) 

0 = 0.036  xR  /b 

X 

(21) 

On  substituting  in  (l8) 


x/pU2  = 0.0288  Rx  1/5  {1  + 1.285  (Ux/U2)}  (22) 

and  on  integrating  the  skin  friction  F of  a width  b of  one  side  of  a 
plate  of  length  1 is  given  by 


Cf  = F/(gpU2lb)  = (2/1)  f (x/pU2)dx 

^ o 

1 • 

= 0.0T2  R-j."  5 (1  + 0.572  (Ul/U2)} 


(23) 
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The  first  term  is  the  steady  state  skin  friction,  while  the  second 
term  is  due  to  acceleration.  The  latter  can  be  written  in  the  form 

F^/;.pUl2b  = M.O'V'i  Rx~  (rh) 

Mince  K contains  U this  form  does  not  lead  us  to  expect  that  the 
'added  mass'  will  be  constant. 

5.  TOWED  PLATE  EXPERIMENT 

It  would  appear  a simple  matter  to  measure  the  resistance  of  a plate 
towed  in  its  plane  in  a fluid  at  accelerating  speed,  and  for  this 
purpose  the  No  2 Ship  Tank  carriage  at  AEW  is  well  suited,  having  a 
reasonably  constant  acceleration  until  approaching  the  ordered  speed. 
In  practice  considerable  care  is  needed  if  the  measurements  are  to  be 
free  of  contamination  from  the  effects  of  roughness,  waviness,  edge 
and  free  surface  effects,  support  interference,  thickness  effect, 
distortion  and  misalignment  of  the  plate  (see  for  example  Hughes 
Reference  13).  In  particular,  if  a thin  plate  is  towed  in  a vertical 
position  by  a support  which  is  clear  of  the  fluid,  that  part  of  the 
leading  edge  away  from  the  support  is  liable  to  deflect  and  possibly 
break.  Dr  Hughes  devised  a special  rig  whicn  enabled  him  to  monitor 
such  deflections,  and  this  is  satisfactory  for  a steady  tow,  but 
greatly  increases  the  towed  mass.  In  the  present  experiment  this  is 
not  acceptable  since,  during  the  accelerated  motion,  the  hydrodynamic 
resistance  will  have  to  be  extracted  from  a gauge  measurement  which 
already  includes  large  inertia  forces. 

This  difficulty  has  been  overcome  by  towing  a plate  curved  into  a 
a U section  (see  Figure  7 and  Plate  1)  the  radius  of  which  is  large 
compared  with  the  boundary  layer  thickness. 

As  can  be  seen  from  section  5-1  the  'Experiment  Details'  this 
configuration  makes  possible  tests  up  to  a Reynolds  Number  of 
6.5  x 106  on  a plate  of  considerable  area  and  rigidity,  yet  having 
a thickness/length  ratio  of  less  than  0.5  per  cent  and  an 
effective  aspect  ratio  (2  x immersed  perimeter/length , as  defined 
by  Dr  Hughes)  of  1.97-  The  radius  of  the  U was  chosen  to  permit  the 
plate  to  enter  the  dock  of  No  2 Ship  Tank , but  there  seems  to  be  no 
reason  to  prevent  plates  of  larger  radius  being  tested  at.  higher 
speeds . 

The  plate  is  a sandwich.  Two  4 mm  sheets  of  plywood,  curved  over  a 
forming  mould  and  glued  together,  are  sandwiched  between  sheets  of 
'Formica'.  This  form  of  construction  ensures  a smooth  finish  free 
from  waviness,  and  eliminates  changes  of  mass  due  to  water  absorbtion. 
With  the  bracing  members  fitted  the  prate  is  extremely  rigid.  By 
limiting  the  length  to  2.438  m (8  feet)  transverse  joins  in  the 
Formica  are  avoided. 

With  the  U configuration  care  is  necessary  to  ensure  vertical 
alignment  as  well  as  fore  and  aft  alignment.  The  latter  was  checked 
during  runs  by  observing  the  water  surface  level  on  the  inside  and 


13 


MO 


outside  of  the  plate.  It  io  also  necessary  to  align  the  swinging 
links  of  the  suspension  to  te  accurately  vertical,  and  the  towing 
link  to  be  accurately  horizontal.  A dynamometer  of  high  stiffness 
was  chosen  so  as  to  ensure  that  the  swinging  links  rer.ai • ed  vertical 
during  the  tow.  The  movement  ol  these  links  was  observe  and 
confirmed  to  bo  negligible. 


5.1.  Experiment  Details 
Curved  Plate 


Length 

Radius  (outer  skin) 

Immersed  perimeter  ( outer  1 
Wetted  aiva  (both  sides) 

Thickness 

Thickness/length  ratio 

2 x wetted  perimeter 
length 


Aspect  ratio 


2.108  m 
711  mm 
2.k05  m 
11 .635  ir2 
1 1 . 5 mm 
0.0  per  cent 

1.97 


Towed  mass 
Surface  fin 
Turbulence 


sn 

timulation 


83.63  kg 
'Formica' 
Sand  strip 
1 ,/8  inch  x 
grains 


0.001k  inch 


Instrumentation 


Carriage  speed  - a.  From  feedback  tachometer  on  carriage 
wheels  recorded  on  U/V  recorder. 

b.  At  steady  speed,  by  time  between  fixed 
points,  electronic  measurement  displayed 
digitally. 

Carriage  acceleration  - from  accelerometer,  recorder  on 

IT  /*  F •» 

v/ v iecot'uer. 

Towing  force  - by  strain  gauged  post  dynamometer,  recorded  on 
U/V  recorder. 

Schedule  of  Runs 


Runs  were  made  at  various  accelerations,  up  to  0.05  g maximum. 

Each  run  consisted  of  an  acceleration  from  rest,  the  acceleration 
being  reasonably  constant,  up  to  various  speeds  (maximum  3.0k8  m/s) 
at  which  the  speed  was  held  constant  for  several  lengths.  Maximum 
Reynolds  Number  attained  6.5  x 106. 

5.2.  Steady  State  Skin  Friction 

The  steady  state  skin  friction  coefficient  is  plotted  in  Figure  8. 
These  values  have  been  corrected  as  in  Reference  13  for  form  drag 
and  turbulence  stimulation  as  follows. 
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Form  drag  correction  (per  cent)  = 3 x thickness/length  ratio 

(per  cent) 

Turbulence  stimulation  correction  = O.OU 

Air  drag  was  measured  by  towing  a duplicate  rig  less  the  in-water  par' 
The  results  are  consistent  with  those  presented  in  Reference  13. 

5.3.  Effect  of  Acceleration 

The  steady  state  drag  (uncorrected)  corresponding  to  the 
instantaneous  speed  has  been  subtracted  from  the  measured  resistance, 
as  has  also  the  inertia  resistance  of  the  plate  and  mountings 
leaving  the  hydrodynamic  drag  due  to  acceleration.  Points  from  a 
number  of  runs  are  plotted  in  non-dimensional  form  in  Figure  <>. 

There  is  a good  deal  of  scatter  on  the  results,  but  it  is  nevertheless 
clear  that  the  momentum  theory  overestimates  the  effect.  At  the  higher 
speeds  the  effect  appears  to  be  negative,  suggesting  that  turbulent 
flow  may  be  suppressed  during  the  accelerated  motion.  However 
this  conclusion  must  be  treated  with  reserve  s'.nce  a discrepancy  has 
been  detected  in  the  measurement  of  the  acceleration.  It  is  intended 
to  reanalyse  some  of  the  points,  but  this  is  unlikely  to  upset  the 
general  conclusion  that  the  momentum  theory  overestimates  the  skin 
friction.  We  may  further  conclude  that  such  increase  in  skin  friction 
as  occurs  is  not  sufficient  to  account  for  the  increase  in  drag  of  the 
submarine.  Thus  attention  should  also  be  directed  to  the  effect  of 
acceleration  on  form  drug. 

6.  CONCLUSIONS 

6.1.  The  increase  of  skin  friction  of  the  laminar  boundary  layer 
a plate  due  to  acceleration  is  well  represented  by 

a.  Blasius  solution  during  the  early  stages  of  the  motion 
ie  before  the  plate  has  travelled  one  length. 

b.  Momentum  theory  during  the  later  stages  ie  after  the  plate 
has  travelled  one  length. 

6.2.  Applying  momentum  theory  based  on  the  seventh  power  law  to  the 
turbulent  boundary  layer  of  an  accelerating  plate  gives  skin  friction 
values  greater  than  those  measured  experimentally.  It  is  possible  that 
acceleration  suppresses  turbulent  flow  in  the  boundary  layer,  but 

some  reanalysis  of  data  is  needed  to  resolve  a discrepancy  in  the 
measurement  of  acceleration. 

6.3.  Skin  friction  increases  do  not  appear  to  be  large  enough  to 
account  for  the  increased  drag  of  a submarine  during  acceleration. 
Attention  should  now  be  directed  to  the  effects  on  form  drag. 
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Appendix  I 


LAMINAR  BOUNDARY  LAYER  COMPUTAT I ON 


Equations  (2),  (3)  and  (M  are  computed,  viz 


at  + 


JU 

3x 


3u 

3y 


U + v 


3 u 

Sy7 


for  the  conditions: 

At  the  surface  of  the  plate  u = v = 0 

Well  away  from  the  plate  u = U,  the  free  stream  velocity 

relative  to  the  plate. 


y 

25 


FREE  STREAM 


i+1 

I 

1-1 


• • • PLATE  SURFACE 

1 1 2 

LEADING 

EDGE 

i-l 

1 

i + l 

25 

TRAILING 

EDGE 

The  value  of  u is  stored  for  each  of  25  values  of  y at  each  of 

25  stations  of  x.  The  x sections  divide  the  length  of  the  plate  into 

25  equal  intervals,  DX  = 1/25,  so  that  the  x co-ordinate  of  the 
ith  station  is  x.  = i DX  where  i = 1 ,25.  (A  station  is  not  required 

at  the  leading  edge,  since  this  is  a singular  point  at  which  no 
computing  is  possible  or  required.)  The  y values  are  equally 

spaced,  DY,  so  that  y.  = (j  - l)DY  where  j = 1,25.  Thus  j = 1 is 

J 
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the  surface  of  the  plate.  The  value  of  DY  may  be  chosen  as  tnought 
appropriate  - in  tr.i^  instance  DY  = 0.001  DX. 

The  solution  is  iterated  at  time  interval  DT  which  may  be  chosen  or 
programmed  as  appropriate  - in  this  instance  DT  = 0.01  seconds. 
f.'Jote  that  too  long  an  iteration  interval  causes  instability  in 
the  computation.) 

In  addition  to  the  main  file  which  stores  all  625  values  of  u, 
four  working  arrays  U1  , U2,  U3,  UL  of  dimension  25  each  are  used. 

The  equations  of  motion  are  evaluated  at  each  x station  in  turn. 

When  the  ith  station  is  being  evaluated  the  existing  velocity  profile-.; 
at  x^  _ 1 , x^,  x^  + ^ are  stored  in  U1 , U2,  U3  respectively.  The 

updated  velocity  profile  at  i is  stored  in  U!t  whence  it  is  written 
into  the  main  file.  Having  completed  the  evaluation  of  the  ith 
station  the  contents  of  U2  are  transferred  to  U1,  those  of  U3  to  U2, 
and  the  velocity  profile  at  i + 2 is  read  from  the  main  file  to  U3 
in  readiness  for  the  evaluation  of  station  i + 1. 


For  eacn  value  of  y at  a given  x station  the  value  of  rr  is  computed 


and  v. 


The  partial  differentials  are  computed  by  simple  differences  for 
values  of  j from  2 to  2k.  (Values  at  j = 1 and  25  are  not  needed 
for  the  computation  of  u since  u = 0 at  j = 1,  and  has  the  free 
stream  value  at  j = 25  - ) 


— = l (U3-  - U1.J/DX  = UX.  sc 
^ X J J J 

!y  ■ i (U2j  - , - U2o  - 2/1>Y 


pony  mil&RIE  TO  DOG  GOES  NOT 
PERMIT  FULLY  LEGIBLE  PROQUCTIOH- 


'Wr 


= (U2.  - 2 U2.  + U2.  J/(DY)‘ 

« J * ■ tJ  il 


The  singularity  at  the  leading  edge  makes  it  impossible  to  compute 
the  velocity  profile  at  the  first  x station.  It  is  assumed  that  the 
boundary  layer  at  this  station  develops  very  quickly  so  that  the 
velocity  profile  is  that  corresponding  to  steady  state  at  the 
instantaneous  speed.  A similar  difficulty  arises  at  station  25 

(the  trailing  edge)  since  it  is  not  possible  to  compute  at  this 

station. 


Strict  precision  required  the  computation  to  be  continued  downstream 

of  the  trailing  edge,  however  is  small  near  the  trailing  edge 

3 li 

and  it  is  sufficiently  accurate  to  base  7^  there  on  the  difference  in 
u between  stations  2k  and  25* 


21 


v is  computed  from  the  equation  of  continuity  (3)  by  integration  of 
using  a trapezium  .rule  and  the  condition  that  at  the  surface  of 

the  plate  v = ~ = 0 

j 

v.  = l }(IK.  1 + UX.  )Di 

J j = 2 J - 1 J 

Since  u is  zero  at  the  surface  of  the  plate  the  skin  friction  at 
station  i is  given  from  (U)  by 

Ti  - “'“j  = 2>i/“ 

Values  of  are  stored  in  an  array  UYO  of  dimension  25,  the  contents 
of  vhich  are  printed  when  all  25  x station  have  been  updated. 

After  printing  output  as  required,  the  time,  distance  travelled 
by  the  plate  velocity  are  updated  in  preparation  for  the  next 
iteration.  It  is  not  necessary  that  the  plate  acceleration  should 
be  constant  - in  this  instance  the  plate  acceleration  is  set  t'' 
zero  when  the  plate  has  travelled  through  three  lengths. 
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Abstract 


Both  full  scale  trial  and  towing  tank  experiment  show  that  the 
'virtual  mass'  in  the  fore  and  aft  direction  is  greater  than 
that  predicted  by  potential  flow  theory.  This  report  examines 
the  increase  in  skin  friction  due  to  acceleration.  A later 
report  will  examine  the  effect  of  acceleration  on  form  drag. 


